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Abstract 



The normal modes of a commensurate monolayer solid may be damped by 
mixing with elastic waves of the substrate. This was shown by B. Hall et 
al., Phys. Rev. B 32, 4932 (1985), for perpendicular adsorbate vibrations 
in the presence of an isotropic elastic medium. That work is generalized 
with an elastic continuum theory of the response of modes of either paral- 
lel or perpendicular polarization for a spherical adsorbate on a hexagonal 
substrate. The results are applied to the discussion of computer simulations 
and inelastic atomic scattering experiments for adsorbates on graphite. The 
extreme anisotropy of the elastic behavior of the graphite leads to quite dif- 
ferent wave vector dependence of the damping for modes polarized perpen- 
dicular and parallel to the substrate. A phenomenological extension of the 
elasticity theory of the graphite to include bond-bending energies improves 
the description of substrate modes with strong anomalous dispersion and en- 
ables a semi-quantitative account of observed avoided crossings of the adlayer 
perpendicular vibration mode and the substrate Rayleigh mode. 

PACS numbers: 68.45.Kg, 68.35.Ja, 63.20.-e 
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I. INTRODUCTION 



In the course of molecular dynamics calculations of the frequency spectrum of com- 
mensurate monolayer solids of nitrogen adsorbed on graphite 0], we became aware of 
paradoxical phenomena in what was expected to be the most ideal and simple regime. At 
low temperatures the center-of-mass one-phonon approximation to the intermediate scatter- 
ing function has a nearly pure sinusoidal oscillation over periods of 20 to 50 ps; however the 
amplitude generally is quite different from that anticipated from equipartition theory for 
harmonic oscillators and from a sum-rule of Hansen and Klein ||. There is some indication 
in the simulation data that the mean-square oscillation amplitude, averaged over hundreds 
of picoseconds, is on the scale expected for oscillator coordinates. 

The only suggestions of such phenomena which we have found in the literature are a com- 
ment by Hansen and Klein H that their sum-rule was satisfied to 10% except at small wave 
numbers where spectral peaks were quite sharp and a comment by Shrimpton and Steele || 
that simulation times much longer than 400 ps would be needed to ensure thermal equili- 
bration of the long-lived Brillouin- zone-center phonons of commensurate Krypton/graphite. 
A nanosecond time scale is inferred from experiments and modeling of the sliding friction 
of incommensurate inert gas monolayers on metal surfaces ||, but the relative importance 
of processes determining such long lifetimes is in dispute ||. We have used perturbation 
theory for the effect of cubic and quartic anharmonicity in the adatom-substrate interaction 
on the lifetime of the adlayer phonons[[7|, [|. With parameters appropriate to commensurate 
Krypton/graphite and to a spherical molecule version of commensurate Nitrogen/graphite, 
the estimated zone-center phonon lifetimes again are on the scale of nanoseconds. Such 
processes are not likely to be the dominant ones in determining the lifetimes, as has been 
appreciated by Mills and his co-workers || 10]. 

The principal mechanism which sets the lifetimes of zone-center phonons in a commen- 
surate monolayer solid is the radiative damping arising because the adlayer normal mode is 
actually a surface resonance that overlaps a continuum of substrate normal modes ||. In 
previous modeling of this process, the substrate was treated as an isotropic elastic medium. 
There was good success in explaining the phonon line-widths (for motions primarily polar- 
ized perpendicular to the substrate surface) of inert gases adsorbed on metals fLOf . The time 
scale for such damping is of the order of picoseconds and the damping is expected to be 
larger for a low-density substrate such as graphite than for high-density metal substrates. 
Inelastic helium atomic scattering experiments for xenon adsorbed on graphite Jl ij and for 
commensurate and incommensurate krypton monolayers on graphite [jl2[ give evidence for 
a strong mixing of the adlayer perpendicular vibration with substrate modes over a range 
of wave vectors and for strong damping of the normal modes at small wavenumbers. 

Here, we develop the elastic substrate theory for the case of adsorption on the (0001) 
(c-axis) surface of a hexagonal substrate. This nominally includes the case of the basal 
plane surface of graphite. However, to mimic the strong anomalous dispersion of the TA^ 
branch of the graphite spectrum [O, HM, the continuum approximation of Komatsu |T5[ and 



Yoshimori and Kitano [16j for the bond-bending energies is adopted. 

There is a nonzero frequency at small wave numbers for perpendicular motions of an 
incommensurate monolayer, mainly determined by the curvature of the adsorbate-substrate 
potential well, and such modes experience both hybridization with substrate modes and 
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damping []9,|17|. For a commensurate monolayer solid, there is a Brillouin zone-center gap 
for motions both parallel and perpendicular to the substrate surface; the radiative damping 
mechanism acts for both polarizations. The long wavelength motions of graphite parallel to 
the surface plane are governed by isotropic elasticity theory, a simplification relative to the 
situation for the (111) surface of face-centered-cubic metals. The large elastic anisotropy of 
graphite between motions parallel and perpendicular to the c-axis has the consequence that 
the strong radiative damping is confined to a much smaller fraction of the adlayer Brillouin 
zone for the branch with parallel polarization than for that with perpendicular polarization. 

The organization of this paper is: Section II describes the models of the interactions 
and the intrinsic dynamics of the decoupled adlayer and substrate. Section III contains the 
formulation of the coupled adlayer and substrate dynamics. Section IV contains the formal 
solution for the adlayer response functions. Some special cases are treated in Section V and 
the results for commensurate monolayers on graphite are presented in Section VI. Section 
VII contains concluding remarks. A summary of our experience with Molecular Dynamics 
simulations for the zone-center modes of the nearly harmonic solid on a static substrate is 
contained in Appendix A. 



II. INTERACTION MODEL AND INTRINSIC DYNAMICS 

The required components are models for the substrate dynamics, the adsorbate- 
adsorbate interactions, and the adsorbate-substrate coupling. The substrate dynamics are 
modeled with an elastic continuum approximation which enables a quite explicit treatment 
at small wave vectors. Apart from an approximation for bond-bending energy terms in the 
graphite substrate |15l|L6| , the dispersion of the substrate normal modes is omitted. The 
adsorbate-adsorbate interactions are taken to be central pair potentials. For the small wave 
number modes, near the Brillouin zone-center gap, the form of the pair potential is not 
crucial to the treatment. Finally, Steele's Fourier decomposition of the adatom-substrate 
interaction is used to make a simple parameterization of the adsorbate-substrate coupling. 
As discussed in Sec. III. D, these choices affect the calculated wave vector dependence of the 
dispersion and damping of the vibrational spectra, but we believe the qualitative features 
of the results are reliable. 

In this paper, the z— axis of the Cartesian coordinate system is taken parallel to the 
c— axis of the hexagonal substrate, the equilibrium configuration of the adsorbed monolayer 
is in the x — y plane, and the wave vector q of the adlayer normal modes is a 2D vector 
lying in this x — y plane. 



A. Intrinsic substrate dynamics 

The hexagonal solid substrate is approximated as an elastic continuum with 5 indepen- 
dent elastic constants. The equations of motion for substrate displacements with components 
Uj are, with the summation convention for repeated indices, 

d 2 ui 

RUi = Ci k l m - . (2.1) 

OXkOXm 
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With hexagonal symmetry, there are only a few nonzero elements in the fourth rank tensor 
Ciki m . In the Voigt notation [18], the subscripts for Cartesian axes are denoted 1 = xx, 2 
= yy, 3 = zz, 4 = yz, 5 = xz, and 6 = xy. Then the independent elastic constants are 
Cn = C22, C33, C44 = C55, C m , an d Ci 3 = C 23 . A further relation derived from the planar 
isotropy is C\i = C\\ — 2Cqq. 

For graphite, the mass density is p = 2.267 gm/cm 3 and the elastic constants are, all 
in 10 11 dyn/cm 2 , C n = 106, C33 = 3.65, C 13 = 1.50, C 44 = 0.40, and C 66 = 44, from 
reference ||19|| . Values for C44 in the recent literature range from 0.325 to 0.47, derived 
from measurements on highly oriented pyrolytic graphite with Brillouin scattering JJIJ an d 
inelastic neutron scattering 
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respectively. We adopt the value C 44 = 0.40 used in 
since it is at the middle of the range and this choice facilitates comparison with lattice 
dynamical calculations [21|J22| . The speed of the Rayleigh wave [pi is nearly equal to that 
of the TAj_ mode for wavevector in the x — y plane and polarization parallel to the c-axis, 
C44/P, because of the large elastic anisotropy of the graphite. 
For the graphite substrate, it is essential to include the strong anomalous dispersion 
13|,14] of the TAj_ mode. This can be accomplished with a continuum approximation to the 
bond-bending energy [llB,|nj by adding the following term to the total substrate energy: 



AE = (A/2) J d 3 r[Vju z + d z V 2 ■ u] 



(2.2) 



where the subscript '2' on the gradient denotes the x — y components. Eq.( [2.2| ) has been 
constructed so that (i) the associated stress tensor is symmetric and (ii) isotropy in the 
graphite basal plane is retained. Then for displacements with spatial dependence 



u(R,t) oc exp(«q ■ R) 



(2.3) 



and wavevectors q in the x — y plane, the generalized elasticity theory for motions in the 
sagittal plane (SP) defined by z and q has the replacement 



C 44 -> C 44 (eff) = C u + Xq 2 



(2.4) 



For motions with shear horizontal polarization (SH), the 'bare' C*4 4 is retained and is denoted 



C44 in this paper. We set A = pK 2 with K = 5.04 x 10 3 cm 2 /sec fitted to the curvature 



of the TAj_ branch observed by inelastic neutron scattering experiments W% . 



B. Adatom-adatom interaction and intrinsic adlayer dynamics 

The adlayer consists of atoms or "spherical molecules" of mass m in a 2D Bravais lattice; 
for commensurate monolayers of krypton or xenon on graphite, it is a triangular lattice. We 
assume that the adatoms interact via a central potential ip and denote by r the projections of 
the equilibrium positions R onto the x — y plane. The analysis of this subsection is for modes 
polarized in the x — y plane; the dominant interaction for the out-of-plane motions arises 
from the adatom-substrate potential treated in Sec. II. C. The normal modes with wavevector 
q have the form 

w j a = w (q) exp(z[q • r ja - w(q)t]) . (2.5) 
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The amplitude and angular frequency are obtained from the solutions of the eigenvalue 
problem 

mw(q) 2 w(q) = D(q)-w(q) (2.6) 
for the dynamical matrix D(q) defined by 

D(q)= E VW[l-cos(q-r,J]. (2.7) 

The small |q| solutions are transverse and longitudinal acoustic waves with frequency pro- 
portional to q. 

The commensurate monolayer on a static substrate has in-plane motions with angular 
frequencies given by 

fi(q) 2 = ^o||+^(q) 2 , (2.8) 

where u>o\\ is the zone-center gap defined in Sec. (2.3). Eq. ( |2.8|) provides the basis for the 
remark at the beginning of this section that the analysis near the Brillouin zone center only 
depends weakly on the form of the pair potential ip. Thus, we adopt the primitive Lennard- 
Jones (12,6) interaction for ip; parameters for krypton and xenon are listed in Table I. 



C. Adatom-substrate interaction 

In Sec. III. A, we assume that the interaction of the adlayer atoms or "spherical" molecules 
j a with substrate atoms J s is given by a sum of central pair potentials: 

$ as = 5X|R, a -R Js |). (2.9) 

ja,Js 

For the case of a static substrate lattice with planar surface, $ as may be transformed fol- 



lowing Steele p4| : 



Astatic = £[K(*;J + E V&ia) eX P(*S • r iJl • (2-10) 

ja g 

The notation in Eq. (|2.10| ) follows that of Sec. II. B: the z-axis is perpendicular to the surface 
and Yj a is the component of R Ja parallel to the surface. The g are the 2D reciprocal lattice 
vectors of the substrate surface. This representation of the static interaction is more general 
than Eq.( |2.9|) since it may include effects of noncentral forces and many-body forces. For 
the following, we truncate the g-sum at the first shell of reciprocal lattice vectors. 

In static substrate models of a commensurate adlayer such as Krypton/graphite with one 
atom of mass m per unit cell, the zone-center modes polarized perpendicular and parallel 
to the surface have angular frequencies given by |25 

<. = ^M^[v (z) + w g (z)]\ z=Zeq 

u 2 0]l =-3g 2 V g (z eg )/m. (2.11) 

where the equilibrium overlayer height is denoted z eq . We use these frequencies to param- 
eterize the dynamic coupling of the adlayer to the substrate, with a further assumption 
which is made explicit in Sec. III. Some values of the frequencies, based on a combination of 
experimental data and modeling, are listed in Table I. 
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III. DYNAMIC COUPLING OF ADLAYER AND SUBSTRATE 



In principle, one might solve the coupled dynamics of the adlayer and the substrate using 
atomistic interaction models for all the constituents. Such calculations have been performed 



21| , |22[| for the normal modes of coupled inert gas-graphite slab systems, but not for the 
effective damping in an adlayer response function. Further, there is only limited knowledge 
of the adatom-substrate corrugation energy. Therefore, we develop a formalism sufficiently 
detailed to show the damping phenomenon and yet one in which the substrate dynamics 
and the adatom-substrate coupling are treated with a few empirical parameters. 



A. Parametric forms 

We must first examine the relation between the descriptions in Sees. II. A and II. C. In the 
former the substrate was treated as a continuous medium with a displacement function u(r, 
z, t); in the latter the atomic discreteness of the substrate was basic to the lateral periodicity 
of the adatom-substrate potential, the amplitudes V g in Eq. (|2.10| ). The formulation for the 
dynamic coupling of the adlayer and the continuum substrate requires a specification of 
where the stress from the adlayer is applied in the substrate. We follow Hall et al. |9j and 
assume it to be concentrated on the surface layer of substrate atoms at height z , displaced 
slightly inward from the boundary z = of the elastic continuum. In the final results, z 
is taken to be vanishingly small. However, an initial distinction between z = and z = z 
is made to bypass complications of 5-function stresses applied precisely at the edge of the 
continuum. 

We specialize immediately to the case where the oscillatory displacements of the substrate 
and the adlayer are represented by 

u Js = u(q, z) exp(>[q ■ r Js - cot)) 

w > = w (q) exp(z[q ■ r ja - cut}) . (3.1) 

The corresponding interaction energy is derived from the second order Taylor series expan- 
sion for the potential $ as of Eq. (|2.9|) 

A$ as = i £ VV0 : (w ja - u Js ) (w ja - u Js y 

= y H VV0 jaJs : [w(q) - exp(zq • [r Js - r ja ])u(q, z)\ 

J s 

x [w(q) - exp(«q ■ [r Js - r jtt ])u(q, z)}* , (3.2) 

where N is the total number of adatoms, the J-sum is assumed to be for atoms in the Zq- 
layer, and Umklapp processes involving reciprocal lattice vectors of the adlayer are neglected. 
Although conclusions about the dispersion based on an atom-atom model for $ as have 
limited generality, such a model was usedj^l],^] for commensurate inert gas / graphite cases 
treated in Sec. VI, so that it is useful to specify the differences in the approach used here. 

Eq.( |3.2| ) may be reduced using a tensor generalization of the analysis which gives 
Eq.( |2.10D . However, in view of the several approximations already made which limit the 
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quantitative accuracy with which the dispersion may be treated, we make one further sim- 
plification and drop the phase factor [27] so that 



1\ r 

A$ QS « - K : J dz 5{z - zo) [w(q) - u(q, z)\ [w(q) - u(q, z)]* , (3.3) 

with the tensor coupling constant K given in dyadic form by 

K = m{ul L zz + c^on [xx + yy}) . (3.4) 

Eq. (|3.4j ) is a parameterized representation for the effect of the adatom-substrate inter- 
action in terms of the Brillouin zone-center gap frequencies discussed in Sec. II. C. It may 
be used to represent the coupling for cases where <3? as is not determined as a sum of pair 
potentials, such as commensurate layers on metals. It is also a way to bypass the incomplete 
understanding of the origin of realistic corrugation amplitudes V g for adsorption on graphite. 
Finally, Eq.( |3.4Q enables a technical simplification in the calculation. When combined with 
the planar elastic isotropy of the hexagonal surface, the problem of coupled adlayer and 
substrate separates into SP and SH motions. 

B. Equations of motion 

The equation of motion for the adlayer normal coordinate becomes 

m uj 2 w(q) = D(q) ■ w(q) + K ■ [w(q) - u(q, z )) ■ (3.5) 

The differential equations for the components of the substrate amplitude u(q,z ) are, for q 
parallel to the x-axis and A tot /N equal to the area per adatom in the commensurate adlayer, 

pu 2 u x {q) = (Cuq 2 - C 44 d 2 )u x (q) 

- iq(C 13 + C A4 )d z u z (q) 

+ (N/A tot ) 8{z - z )K 0xx [u x (q) - w x (q)} 
pu 2 u z {q) = (C u q 2 - C 33 d 2 )u z (q) 

- iq(C 13 + C A4 )d z u x (q) 

+ (N/A tot ) 5(z - z )K 0zz [u z (q) - w z (q)} . (3.6) 

and 

puj 2 u y (q) = (C m q 2 - C ( $d 2 z )u y (q) 

+ {N/Atot) 5(z - z )K 0yy [u y (q) - w y {q)} (3.7) 



C. Boundary conditions 

As the boundary conditions on the substrate displacement function u(q, z), we take 
the z = boundary to be a free surface where the following components of the stress tensor 
vanish: 
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T z p\z=o = 0, (3 = x,y,z. (3.8) 
The components of the stress tensor are given by 

I' a, /3 — c a/3kl d k ui , (3.9) 



using the 4- index form of the elastic constants. For the displacement function of Eq. (|3.1| ), 
with q parallel to the x-axis and returning to Voigt notation, Eqs.(|3.8|) become (with z = 0): 

d z u x (q,z) + iqu z (q,z) = 
d z u y (q, z) = 

d z u z (q, z) + iq(C 13 /C 33 )u x (q, z) = . (3.10) 

Finally, the theory of the adlayer response involves substrate motions driven by an initial 
displacement of adlayer atoms. Then, deep in the substrate, z — > — oo, the disturbance 
created by the adlayer must decay exponentially or take the form of an 'outgoing' wave. 
This becomes a requirement that the solutions of Eqs.( |3.6| ) and ( |3.7| ) for z < zq have the 
form exp(— K\z\) or exp(zK\z\) with K > - see Sec.IV.B. 



D. Comments 



Eqs.( |3.5| ) to ( |3.7| ) generalize the treatment of Hall et al. in two ways: the substrate 
is an anisotropic elastic continuum and there are driving terms which arise from the coher- 
ent addition of lateral force terms for the commensurate adlayer. There is an increase in 
complexity beyond their treatment, but, as shown in Sec.V.C, quite simple results are again 
obtained at the Brillouin zone center. 

We summarize the approximations that have been made which have serious consequences 
for the treatment of the dependence of the mode damping on the wave vector: 



1. A distinction is made between the edge of the elastic continuum at z = and the 
height Zq where the adlayer stress is applied. However the limit Zq — > is taken in the 
analysis. 

2. The anomalous dispersion of the TAj_ branch of the graphite substrate is approximated 
in the elastic continuum description by using the effective elastic constant G^eff) 
defined in Eq.(|2.4j). This replaces C44 in Eqs.(|3.6|). It leads to a large shift in the wave 
number where the Rayleigh wave of the graphite hybridizes with the u± adlayer mode 
and improves the agreement with the HAS experiments |lljl2 . 



3. If the dynamical matrix D(q) is dropped from the adlayer equation of motion, there 
is only a small effect for small wave numbers. 



The approximation in dropping certain phase factors to obtain Eq.( [3.3|) is accurate at 
small wave numbers; however it omits a q— dependence of the dynamic coupling of the 
adlayer and substrate p7 . 
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IV. CORRELATION FUNCTIONS 



The response of the adlayer in the presence of the substrate is characterized using the 
time Fourier transform of correlation functions of displacement amplitudes defined by 

S aa (q,t) = (W a (q,t)W a (q,0)), (4.1) 

where a — x,y, z. The initial conditions on the displacements are 

W a (q,t = 0) = W a0 

W a (q,t = 0) = 0, (4.2) 

with zero for t < 0, and the substrate is initially unperturbed and static. Then the Fourier 
transform for Eq. (|3.5| ) is generalized to 



CO 







exp(iu)t)W a (q, t) dt = —u 2 w a (q, uj) + iu)W a Q , (4.3) 



using 



1 r°° 

W a (q,t) = — w a (q,uj)exp(-iujt)duj . (4.4) 

Z7T J-oo 

In this and the following sections, the dependence on wave number q has been omitted from 
the notation, to reduce the complexity of the formulae. 



A. Green's function solution 

The solution to the set of coupled dynamical equations posed in Sees. III. B and III.C is 
conveniently stated in terms of the values at z = z' = z of a set of Green's functions g Ql g(q, 



a>, z \z') satisfying the following set of equations |£J 

(L • g) ajS = 5 a(3 5(z - z) , (4.5) 

where the 3x3 differential tensor L has the following nonzero elements 

L xx = pu 2 - C u q 2 + Cud 2 
L zz = puj 2 - C M q 2 + C 33 d 2 z 
L yy = pu 2 — C§%q 2 + Cud 2 

L zx = L xz = tq(C 13 + C u )d z . (4.6) 



The boundary conditions at z = based on Eqs.( |3.1(3D are, for a = x,y, z, 



d z g xa + iqg za = 

dzQya = 

c ■ 

d z g za + iq-p^-Qxa = . (4.7) 

^33 
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The problem separates so that the functions g xy , g yz , g yx , and g zy vanish . 
Then with the definitions 

X x = (N/A tat )rnuj^ 

X z = (N/A tot )muj 2 L , (4.8) 

the functions Uq,(z ) are given in terms of the functions g a p = g a /3(zo\zo) by 

u x = g xx X x (u x - w x ) + g xz X z (u z - w z ) 
u z = g zx X x (u x - w x ) + g zz X z {u z - w z ) 

u y = g yy X x (u y - w y ) . (4.9) 

A formal solution for the driving terms in the adlayer equations of motion is 

u x -w x = b n w x + b u w z 
u z -w z = b 2 iw x + b 22 w z 

Uy-Wy = Wy/ ^yy X X ~ l) , (4.K)) 

where the coefficients are given by 

feu = \g Z zX z - iyw h 

bn = -9xzX z /W b 
hi = -g Z xX x /W b 

b 2 2 = [g xx X x -l]/W b (4.11) 

and 

W b = [g xx X x - l][g zz X z - 1] - g xz X z g zx X x . (4.12) 

The adlayer equations of motion then become 

[u 2 - iul(q) + ul\\b n \w x + ^011^12^2 = iuW x0 

UQj_b 21 w x + [uj 2 + uj 2 j_b 22 ]w z = iu}W z0 (4.13) 

and 



2 u?(q) + (u 2 J[g yy X x - !})}w y = iujW y0 . (4.14) 



The u)e and u>t are the frequencies of longitudinal and transverse polarization, respectively, 
in the intrinsic adlayer dynamics, Eq.( [2.6|) , and the x— axis is taken to be a high symmetry 
direction, TM or TK, of the adlayer Brillouin zone. 
The solutions are used to form 

S aa (q,u) = \w a (u,q)\ 2 (4.15) 

with initial condition 

W fiQ = 6 aP . (4.16) 

The damping of the adlayer normal modes manifests itself as broadened peaks in S(q, uS) 
as a function of u at fixed q. Generally, peaks in S(q, uj) may be assigned as derived from 
the intrinsic adlayer frequencies or from the Rayleigh wave of the substrate. As shown in 
Sec. IV. B, the radiative damping mechanism operates for sufficiently small q. 
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B. Evaluation of the Green's functions 



The problem separates into analysis of the SP and SH motions with the coupled x — z 
equations and ^/-equation, respectively. 



1. Sagittal plane 



The solution is very similar to one given by Dobrzynski and Maradudin |23[ for the 
Green's function of a hexagonal elastic half space. We seek solutions of the homogeneous 
versions of Eqs. (|4.5|) which have exponential dependences on z: 

gps ~ exp(az) . (4.17) 



With the definitions 



71 = [pu 2 - C n q 2 ]/Cu 
74 = [pu 2 - C 44 q 2 ]/C 33 
ox = 7i + 74 + g 2 [(Ci 3 + C 44 ) 2 /C 33 C 44 ] 



(J2 = V °l - 4 7i74 , (4.18) 
there are two inverse length scales ay given by: 

a 2 = [- ai +a 2 ]/2 

a 2 2 =[-a l -a 2 }/2. (4.19) 

For z < z', the roots of Eq. (|4.19|) are chosen to give damped or outgoing waves according 
to whether a-j is real or imaginary ||29||. Denote the longitudinal acoustic and transverse 
acoustic modes for wave vector in the x — y plane by LA (SPy) and TA^, respectively, and 
the LA mode for wave vector along the z— axis by LA^. The corresponding speeds in the 
long wavelength limit are 



CLA = 


s/Cn/p 


C T A = 


Jc 44 /p 


CLAz = 


J C 33 /p 



(4.20) 



The choice of roots for Eqs.( |4.19| ) is then 



and 



an = , ^ > c LA q, 

ai = |«i|, uj<c LA q, (4.21) 



a 2 = — i\ot2 \ , ^ > cta q , 

«2 = |a 2 | , ^<c TA q- (4.22) 
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According to Eqs.( [4.5|) , the Green's functions are coupled in pairs (g xx , g zx ) and (g xz , 
g zz ). Then, for Eq.( 4.17f ) we have 



9xp{aj) = ifjg z/ 3(oij) , 
with the proportionality factor fj defined by 

fj = -qaj(C 13 + C^)l\C^{a) + 7l )] . 



The solution of the homogeneous form of Eqs.( [4.5| ) in the range z < z' is 

9x0 = aexp(aiz) + 5exp(a 2 ;z) 

g x /3 = iafx exp(aiz) + ibf 2 exp(ct 2 z) , 

and in the range z 1 < z < is 

g zl3 = [A + exp(a!z) + A„ exj)(-aiz)) 
+ [B + exp(a 2 z) + B_ exp(-a 2 ^)] 

9xp = ifi [A+ exp(a 1 z) - iA_ exp(-aiz)] 
+ if 2 [B + exp(a 2 -2) - iB_ exp(— a 2 z)) . 



(4.23) 
(4.24) 

(4.25) 



(4.26) 



The six coefficients a, b, A + , A_, B + and B_ are determined from six equations: the z = 
boundary condition Eqs.( [4.7|) , the continuity of g x p and g z p at z = z', and the matching of 
the discontinuities in the first derivatives at z = z' to the strengths of the <5-functions. The 
latter equations are 



C 4 4 [d z g, 
C33 [d z g 



xx z- z 



- d z g 

- d z g zz 



XX z = z' 



1 

1 . 



(4.27) 



The z— derivatives of g xz and g zx are continuous at z — z'. Completion of the explicit 
solution for the Green's functions then is an exercise in linear algebra. 

The solutions for g a/3 (z'\z') with z' — > can be given in compact form using the defini- 



tions: 



a u = q + ocifi 

a 12 = q + a 2 f 2 

a 21 =a x - g(Ci3/C 33 )/i 

a 22 = ol 2 - q(C 13 /C 33 )f 2 



and the Wronskian 



W a = a n cL22 - ana 2 i . 
The Green's function components are 

9zx = -«( Q 21 - CL 22 )/W a C AA 

g X x = (M21 - fia 22 )/W a C 4 4 
g xz = -2(/20H - fia 12 )/W a C 33 
g Z z = (012 - a n )/W a C 33 . 



(4.28) 
(4.29) 



(4.30) 
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There are three q— ranges: I, q < u/cla, where both transverse and longitudinal substrate 
waves are involved in the damping; II, uj/cla < q < uj/cta, where only the transverse waves 
are involved; and III, uo/cta < Q, where there is no radiative damping. Characteristic values 
for uj are ui \\ for the parallel polarization mode and u> ± for the perpendicular polarization. 

2. SH mode 

Define 

al = (C 66 q 2 - pu 2 )/C<$ - (4.31) 
The speed of the transverse elastic waves in the x — y plane, denoted the SH mode, is 



c S H = ^C w /p, (4.32) 
and the choice of root of Eq. ( |4.31|) is 

a 3 = -i\a 3 \ , u > c SH q, 

a 3 = |a 3 |, uj<c SH q- (4.33) 



Then the solution for g yy has the formpO 



g yy = Yi exp(a 3 z), z < z 

= Y 2 cosh(a 3 2:), z' < z < , (4.34) 



where the z = boundary condition, Eq. (|4.7|) , has been used. The coefficients Y\ and Y 2 
are obtained from the continuity conditions at z — z' 

g yy (z = z'+) - g m (z = z'-) = 
C*44^ [dz9yy\z=z'+ ~ d z g yy \ z=z '-] = 1 . (4.35) 

The solution for Y 2 is 

Y 2 = -exp(a 3 z')/(a 3 c£>). (4.36) 
Then, with z - > 0, the value of g yy entering in Eq. ([4.14|) is 

g yy = -l/(a 3 C$). (4.37) 

V. SPECIAL CASES 

We discuss here three special cases where the present formalism overlaps with other work. 
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A. Rayleigh wave 



The frequency (speed) of the Rayleigh wave of wave number q is the root of W a = 0, 
for the Wronskian defined in Eq.( |4.29|) . In the limit q — > 0, this reproduces the result of 
Dobrzynski and Maradudin f23f . As noted by others [19], in the small— q limit the speed 



of the Rayleigh wave of the graphite basal plane surface is only 0.02% smaller than cta- 
The solution for the Rayleigh wave frequency at finite q with the effective elastic constant 
Eq.( [2.4| ) is formally the same, but the quantitative results change somewhat. With the 
parameters used here, the Rayleigh frequency is 0.1% smaller than the TAj_ frequency at 
q = 0.3 A -1 and 0.8% smaller at q = 0.6 A -1 . Even so, the difference between the Rayleigh 
frequency and the TAj_ frequency remains much smaller than the 8% difference found for 
the case of a Cauchy isotropic elastic solid. 



B. Isotropic elastic medium 

The formalism reduces to the case treated by Hall et al. M by choosing 



W0[| 


= 


C n 


= c. 


C44 


= c, 


C*13 


= c 



33 
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(5.1) 



and examining the structure of the response function Szz(q,u) for fixed q. Results for 
the peak frequencies and full-widths at half-maximum for the damped peaks of Szz{q,u) 
are shown for a model of Xe/Ag(lll) in Figure 1. For Figure 1, we extended slightly the 



original calculation of Hall et al. P,|3T]| using the parameters luq± = 2.8 meV (0.67g THz) 
and p = 10.635 gm/cm 3 and omitting adatom - adatom interactions (ip = 0). The effective 
elastic constants Cn = 17.7 and = 2.86 (10 11 dyn/cm 2 ) were fitted to the calculated 



speeds of longitudinal and transverse sound for the Ag(lll) surface |32 



Qualitatively || |It]], the peak frequencies follow trajectories characteristic of an avoided 
level crossing of the substrate Rayleigh wave and the u± adlayer mode at q ~ 0.3 A -1 . 
For q < uj±/cta, the u±— mode is damped and there is a sharp resonance at a frequency 
somewhat reduced (the avoided crossing phenomenon) from that of the bare Rayleigh wave. 
At q ~ uj±/cla, near 0.1 A" 1 , there is an additional contribution to the damping and a 
perturbation to the peak frequency derived from u±, a phenomenon that has been termed a 
van Hove anomaly f33fl . The branch which is the Rayleigh mode at small q approaches Uq± 
at large g, but is still 7.5% below that limit at 0.4 A -1 . 

A novel feature occurs for the present choice of parameters: there is only one sharp 
resonance at small q, but at sufficiently large q there are two sharp resonances. One is 
derived from the Rayleigh mode and one from u±. The second sharp resonance arises 
because the upper 'repelled' frequency lies between the bare substrate Rayleigh frequency 
cr q and the continuum of substrate frequencies that begins at cta <?• There is a 6% difference 
between cr and cta in this model. That there is a signature of the substrate Rayleigh wave 
at wave numbers both above and below the avoided crossing has been considered a notable 
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phenomenon in helium atom scattering from adsorbed monolayers ||34j| . We do not find 



the corresponding effect in the calculations for adsorbates on graphite, Sec. VI, apparently 
because there the increment between the Rayleigh frequency and the bulk continuum is 
quite small. 



C. Small q-limit 

In the q — > limit, the results of Sec. IV have simple explicit forms. The coefficients g zx , 
9xz, bi2, and 6 2 i vanish, so that the w x , w y , and w z motions are decoupled. The remaining 
Green's function components become (for u > 0) 

9xx = -l/[Cu<*2] = -i/lpCTA^] 

9zz = -l/[C 33 «i] = -i/[pc LAz ui\ 

9yy = 9xx ■ (5-2) 



Then, defining 



r x = K/[pcta] 

T z = X z /[pCLAz] ■ (5.3) 



the spectral functions are 



2\2 



s «l°rf- ^ + n-<)> + <ri • (5 ' 4) 

and Syy(0,u) = Sxxi^,^)- 

Approximate expressions for the full- widths at half- maximum for Sxx an d Szz, respec- 
tively, are 

5lu z ~ T z . (5.5) 

Eqs.( |5.3j ) show that the damping is enhanced for lower density substrates if the other pa- 
rameters remain similar. This indeed is the trend found in comparing the damping of the 
u±— modes of Xe/Ag(lll) and Xe/graphite, see Sec. VI. Eqs.(|5.5|) are accurate for weak 
damping; the results reported in Sec. VI are obtained with the full formalism of Sec. IV and 
include self-consistent solutions for cases with strong damping. 

Using the A r 2 /graphite parameters in Table I, we find T x /u\\ ~ 0.25 and an estimate of 
3 ps for the decay time. This supports the assertion in Sec. I that the radiative damping 
mechanism is the dominant process determining the lifetime of the zone-center mode. 
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VI. COMMENSURATE MONOLAYERS ON GRAPHITE 



We present applications of the elastic substrate theory of radiative damping to com- 
mensurate monolayers of Xe/graphite and Kr /graphite and also compare to the inert gas / 
graphite slab frequency spectra calculated by DeWette et al. ■ Although the lateral 

interactions are the same as in that work, we have adjusted the frequencies ujq\\ and ujq± to 
incorporate more recent information, so that there are quantitative differences which arise 
from differences in the interaction models. 

Figure 2 shows the results for the Kr/graphite a/3— commensurate monolayer and Figure 
3 shows the results for the corresponding Xe/graphite case. The direction of q is along the 
TK axis of the adlayer Brillouin zone. The ojq±— frequency is marked as a solid horizontal 
line in both graphs and dotted and dot-dash lines denote the thresholds of bulk graphite 
continua based on the SPn and TAj_ modes, respectively. Lateral interactions, with the 
parameters of Table I, are included following the discussion of Sees. II. B and IV.A. The 
plotted points are the derived peak frequencies of Sxx and Szz, as noted, with widths of 
damped peaks indicated by error bars. Cases where an error bar coincides with a substrate 
threshold denote a local minimum for the response function, without a full decrease to half 
the peak height. 



Experiments with Helium Atomic Scattering (HAS) |TT] , 12 1 indicate there is a strong 
damping of the uj±— mode at small q and show a strong perturbation for q « 0.25 — 0.3 
A -1 where the TAj_— mode of the bare graphite crosses ujq±. The extrapolated crossing 
using the initial slope of the TA^— mode is q ~ 0.4 — 0.5 A~ x . However, including the 
strong anomalous dispersion of the TA^— branch with the prescription in Eq.( |2.4j) , leads to 
a semi-quantitative account of the crossing. 

Second, we examine the radiative damping for the in-plane zone-center gap wq\\. The 
region of strong damping for peaks of Sxx(q, w) is confined to region I defined in Sec.IV.B.l, 
i.e., to the left of the substrate SPy threshold shown in the Figures. In region II, between the 
SPy and TAj_ thresholds, the widths of the peaks in Sxx are small but finite and correspond 
to lifetimes on the scale of ns. The large elastic anisotropy of the graphite makes region 
I much smaller than for the isotropic elastic medium: using the parameters for Ag(lll) 
in Sec.V.B the ratio cta/cla is 0.38, but for graphite it is less than 0.1 for q < 0.3 A -1 . 
Another manifestation of the large elastic anisotropy is that the elliptical polarization of the 
Rayleigh wave nearly degenerates to a transverse (z) polarization, with only weak coupling 
to in-plane motions of the adlayer. 

Third, the damping of the uo±_ branch is very strong in both regions I and II. In contrast 
to the model for Xe/Ag(lll), Figure 1, we obtain only one sharp resonance (denoted by +) 
in Szz for region III. The Figures show a large shift of the resonant frequency in region III 
relative to the value for the static substrate used as input to the calculation. 

Finally, we compare the resonant frequencies themselves with the atomistic calculations 
of DeWette and coworkers [21 , 22 1 at large q for a test of the size of effects of the neglected 
spatial dispersion. The most significant discrepancy is that 'deflection' of the trajectories in 
the region of the avoided level crossing is much larger for the continuum calculation than 
in the atomistic calculation. In the elastic continuum theory, the shift remains on the order 
of 10% to the largest q of the calculation; this is a 50% larger shift than in the atomistic 
calculations. The dispersion with q of the peak l u x ' (x and o) of Sxx, from effects of 
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adatom-adatom interactions, is similar to that in the atomistic calculations. The present 
calculations show the uj x branch crossing the TAj_ branch (actually, the Rayleigh mode), as 
do the atomistic calculations. 

The main previous test of the elastic continuum theory of radiative damping was for the 
damping of incommensurate inert gas adlayers on Pt(lll) [Kj, where the formalism tended 
to underestimate the zone-center damping. The factor of ten in the mass density of the 
substrate between graphite and platinum has the effect of making the damping much larger 
for Xe/graphite. This was anticipated by Toennies and Vollmer in their discussion of 



the rather broad peaks for the u±— mode in the HAS inelastic scattering experiments. 



VII. PROSPECTS 

These calculations show that the radiative damping mechanism proposed by Hall et al. 
has a major effect on line-widths which may be observed in inelastic scattering experiments 
from commensurate adlayers on graphite. The understanding of the coupling of the com- 
mensurate layer to the substrate is more advanced for substrates such as graphite than for 
metallic substrates. Thus, adsorbates on graphite may be good subjects for detailed further 
study. It would be of interest to determine whether there are related effects of substrate 
dynamics on the monolayer fluid which would disrupt the long-time tails seen in molecular 
dynamics simulations of the N^/graphite fluid for 1 to 10 ps. 

Another question is how to relate the size of the avoided level crossing of an adlayer mode 
and the substrate Rayleigh mode to adlayer-substrate coupling constants. Comparison of 
our elastic continuum results to the model calculations of DeWette and co-workers indicates 
that there are significant effects of spatial dispersion to be included. This might be explored 



in future work based on a technique such as lattice Green's functions|10[, now that the 
elastic continuum theory is in place. For the damping at intermediate and large wave 
numbers, where the radiative damping mechanism becomes small, a treatment of the more 



conventional anharmonic damping will be needed [10 . 

The large differences in the damping of parallel and perpendicular motions for the com- 
mensurate monolayer on graphite seem well-based and may have ironic consequences. The 
HAS experiments for such monolayers could have more prominent inelastic peaks for the 
parallel than for the perpendicular motions, in spite of the role of polarization considerations 
in the coupling to the helium atom to the adlayer. 
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APPENDIX A: MOLECULAR DYNAMICS TREATMENT OF THE 

ZONE-CENTER MODE 



We summarize the considerations for the center-of-mass one-phonon approximation to 
the intermediate scattering function. Such calculations for a static-substrate model of the 
commensurate ^/graphite monolayer actually were the starting point of the present paper. 

Define the collective coordinate 

1 N 

*) = m £ ex P( zc t ■ R i) <1 ■ 5R ^ ■ ( A1 ) 

The corresponding intermediate scattering function is 

F 1 (q,t) = (X(q,t)X(q,0)), (A2) 

with Fourier transform 

1 r°° 

Si(ti,u) = — Fi(q,t)exp(-«jt)dt. (A3) 

This definition omits a Debye- Waller factor present in the analysis of Maradudin and Fein 
||. The brackets denote a thermal average at temperature T. Hansen and Klein || derived a 
sum-rule for the classical-mechanics limit pertaining to the molecular dynamics calculation: 

/oo 
u 2 S 1 (q, uo) duo = k B Tq 2 /Nm . ( A4) 

-oo 

In the molecular dynamics calculations at a reciprocal lattice vector q = r, corresponding 
to the Brillouin zone-center "0" , there are long time-series where F\ appears to have a purely 
harmonic variation: 

Fx(q,t) ~ A F cos(cu t) ■ (A5) 

Then the sum-rule gives 

Ap = k B Tr 2 /Nmio 2 . (A6) 

Eq. ( |A6| ) is also obtained from the equipartition of energy for an oscillator coordinate Q a . 
Take r to be along an axis a: 

±Mu 2 (Q 2 a ) = h B T. (A7) 

This leads to the same expression for Ap, using M = Nm. The oscillator representation 
has the additional consequence that the fluctuations in Ap may be derived from the fourth 
moment of the Gaussian which leads to Eq. ( |A7| ) : 

(A 2 F )=3(A F ) 2 . (A8) 

Eq. ( |A6|) is badly violated in the results we presented in Fig. (13) of ref.|l[ for two re- 
ciprocal lattice vectors of the herringbone lattice of ^/graphite. At q = 1.703xA _1 , the 
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amplitude is much less than the equipartition result (N = 224) and at q = 2.95yA~ 1 , it is 
much greater. While those results were obtained for time series of ~ 40 ps, extending the 
time series to 400 ps did not improve the agreement with Eq. ([A6|) by much. Results for 
this pair of reciprocal lattice vectors still had no simple relation to each other, as individual 
amplitudes were much less or much larger than expected from equipartition theory. There 
was some indication that Eq.(JA^) was being satisfied. 

The failure of the MD to match to equipartition strength for good harmonic low tem- 
perature solids was observed for the complete monolayer p = 1.0 and for the submonolayer 
p = 0.5. There also were time intervals where the apparent amplitude Ap 'switched' to 
another value; there was a beating pattern in F\ composed of two (and sometimes three) 
closely spaced frequencies. The phenomena were seen both in the center-of-mass and in the 
atomic coordinate versions of the intermediate scattering functions. The time scales were in 
the range 50 — 200 ps. Application of the theory of radiative damping shows that the dy- 
namic substrate has an important role in achieving thermal equilibration of the zone-center 
modes already at much shorter times. 
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APPENDIX: FIGURE CAPTIONS 



1. Model of incommensurate Xe/Ag(lll). The frequencies (in THz) for peaks in 
SzziQ-,^) are plotted as a function of wave number q (in A~ x ). The error flags on 
the symbols A denote the frequencies at the half-maxima for given q. The points + 
denote sharp resonances. The solid horizontal line denotes the 'bare' frequency u>o± 
and the three straight lines with successively smaller slope are the dispersion relations 
for the longitudinal and transverse acoustic modes and the Rayleigh wave at the silver 
surface. The substrate is modeled as an isotropic elastic continuum. See Sec.V.B for 
the parameters used. 

2. Commensurate Krypton/graphite. The frequencies (in THz) for peaks in Szz{q, and 
Sxx(q, as noted, are plotted as a function of wave number q (in A" 1 ) directed along 
the TK axis of the adlayer Brillouin zone. The A and + denote frequencies derived 
from the peaks in Szz, while x, □, and o denote frequencies derived from peaks in Sxx- 
The error flags denote the frequencies at the half-maxima in the response functions. 
The peaks denoted by + and o are sharp resonances. In the g-range where there are 
no error flags on the x, the widths of the peaks are finite but too narrow to be visible 
on the scale of this graph. The values of c^oil and u ± are indicated by dashed and 
solid horizontal lines respectively. The thresholds of the TA^ and SPii continua of the 
graphite are indicated by dot-dash and dotted lines respectively. For the parameters 
of the interaction model see Table I. 

3. Commensurate Xenon/graphite. Identifications as in Fig.(^). Note the crossing of the 
uj x and Rayleigh branches with no apparent deflection. 
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TABLES 



TABLE I. Interaction parameters and frequencies of perpendicular and parallel adlayer motions 
at the Brillouin zone center for a commensurate v3 lattice on a static graphite substrate; frequencies 
are in meV. a 



Case 


e b 


a b 


m b 






Kr /graphite 


159 c 


3.60 c 


139.2 


4.1 d 


1.0| 


Xe/graphite 


228-f 


3.97 f 


218.0 


3.09 


0.57^ 


N2 / graphite 


95 h 


3.7 h 


46.76 


6.0* 


1.6| 



a) The conversion factor to THz is 1 meV = 0.242 THz. 

b) Lennard Jones (12,6) e in Kelvin and a in A; mass in 10~ 24 gm/mol. 

c) parameters from |22| . 



d) from HAS experiment 

e) Value based on model calculation with enhanced corrugation 



f) from pi 



g) from HAS experiment |TT 



h) Values for 'spherical nitrogen' to match 3D critical temperature and nearest neighbor 
spacing in ground state solid. 

i) from inelastic neutron scattering [^6| on orientationally ordered herringbone lattice. 
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